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ABSTRACT 

We use a series of Monte Carlo simulations to investigate the theory of galaxy-galaxy 
lensing by non-spherical dark matter haloes. The simulations include a careful ac- 
counting of the effects of multiple deflections on the galaxy-galaxy lensing signal. In 
a typical observational data set where the mean tangential shear of sources with red- 
shifts ~ 0.6 is measured with respect to the observed symmetry axes of foreground 
galaxies with redshifts zi ~ 0.3, we find that the signature of anisotropic galaxy-galaxy 
lensing differs substantially from the simple expectation that one would have in the 
absence of multiple deflections. In general, the observed ratio of the mean tangential 
shears, 'y~^{0)/j^{d), is strongly suppressed compared to the function that one would 
measure if the intrinsic symmetry axes of the foreground galaxies were known. De- 
pending upon the characteristic masses of the lenses, the observed ratio of the mean 
tangential shears may be consistent with an isotropic signal (despite the fact that the 
lenses are non-spherical), or it may even be reversed from the expected signal (i.e., 
the mean tangential shear for sources close to the observed minor axes of the lenses 
may exceed the mean tangential shear for sources close to the observed major axes of 
the lenses). These effects are caused primarily by the fact that the images of the lens 
galaxies have, themselves, been lensed and therefore the observed symmetry axes of 
the lens galaxies differ from their intrinsic symmetry axes. We show that the effects 
of lensing of the foreground galaxies on the observed function ■y^ {0) / -f~ (6) cannot be 
eliminated simply by the rejection of foreground galaxies with very small image ellip- 
ticities, nor by simply focusing the analysis on sources that are located very close to 
the observed symmetry axes of the foreground galaxies. We conclude that any attempt 
to use a measurement of ■y~^{9)/j~{6) to constrain the shapes of dark matter galaxy 
haloes must include Monte Carlo simulations that take multiple deflections properly 
into account. 
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1 INTRODUCTION 

Galaxy-galaxy lensing is a form of weak gravitational lens- 
ing in which background galaxies are systematically lensed 
by foreground galaxies. Brainerd, Blandford & Small (1996; 
BBS) published the first statistically-significant (4(t) detec- 
tion of this effect using a small data set that consisted of 
439 foreground galaxies, 506 background galaxies, and 3202 
foreground-background galaxy pairs. Since this early work, 
galajcy-galaxy lensing has been detected with high precision 
using various data sets, most of which contain millions of 
foreground-background galaxy pairs. These high-precision 
detections have allowed direct constraints to be placed on 
the nature of the dark matter haloes that surround the lens 
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galaxies, as well as on the bias between mass and light in 
the universe (see, e.g., Fischer et al. 2000; Guzik & Seljak 
2002; Hoekstra, Yee & Gladders 2004; Hoekstra et al. 2005; 
Sheldon et al. 2004; Heymans et al. 2006; Kleinheinrich et 
al. 2006; Mandelbaum et al. 2006ab; Mandelbaum, Seljak 
& Hirata 2008; Limousin et al. 2007; Parker et al. 2007; 
Natarajan et al. 2009; Tian et al. 2009). 

Observations of galaxy-galaxy lensing by field galaxies 
have shown: [1] at fixed luminosity, the haloes of red (early- 
type) galaxies are more massive by a factor of ~ 2 than 
the haloes of blue (late-type) galaxies (e.g., Guzik & Seljak 
2002; Kleinheinrich et al. 2006; Sheldon et al. 2004; Mandel- 
baum et al. 2006a), [2] the haloes of high- luminosity galaxies 
are more massive than the haloes of low-luminosity galaxies 
(e.g., Sheldon et al. 2004; Mandelbaum et al. 2006a), and [3] 
the dark matter profiles of the haloes are consistent with the 
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spherically-averaged Navarro, Frcrik & White (NFW) profile 
(Navarro, Frenk & White 1995, 1996, 1997; e.g., Heymans 
at al. 2006; Hoekstra et al. 2004, 2005; Kleinheinrich et al. 
2006; Mandelbaum et al. 2008). In other words, observa- 
tions of galaxy-galaxy Icnsing by field galaxies have yielded 
a picture of luminous gala^xies and their dark matter haloes 
that is broadly consistent with the expectations of galaxy 
formation in the context of the cold dark matter (CDM) 
model. 

Despite the popularity of the spherically-averaged NFW 
density profile, CDM haloes axe not spherical. Rather, CDM 
haloes are triaxial and the degree of flattening increases with 
halo viral mass (e.g., Warren ct al. 1992; Jing & Suto 2002; 
Bailin & Steinmctz 2005; Kasun & Evrard 2005; AUgood 
et al. 2006). In principle, galaxy-galaxy lensing should be 
able to provide constraints on the shapes of the dark matter 
haloes of field galaxies, since a non-spherical weak lens will 
produce an anisotropic shear pattern. Consider an isolated 
weak galaxy lens with a non-spherical dark matter halo (i.e., 
a halo that, in projection on the sky, has an elliptical surface 
mass density). For fixed source redshift and fixed angular 
distance from the lens, sources that are located closer to 
the major axis of the lens will experience greater shear than 
sources that are located closer to the minor a^xis of the lens. 
If the halo of the lens can be approximated as a singular 
isothermal ellipsoid with projected ellipticity ehaio = 0.3, the 
shear experienced by sources nearest the minor axis of the 
lens will be ^ 80% that of the shear experienced by sources 
nearest the major a^cis of the lens (see, e.g., Brainerd & 
Blandford 2002). Although small, such an anisotropy in the 
galaxy-galaxy lensing signal should be observable provided 
that, in projection on the sky, mass and light arc reasonably 
well aligned within the lens galaxies. 

Weak lensing by galaxy clusters in the Sloan Digital 
Sky Survey (SDSS; e.g., Abazajian et al. 2009 and refer- 
ences therein) has shown that the dark mass associated with 
galaxy clusters is non-spherical and has a projected axis ra- 
tio of b/a = 0.48lo Jg (Evans & Bridle 2009). The detection 
of non-spherical haloes by galaxy-galaxy lensing has, how- 
ever, proven to be more problematical. In a study of galaxy- 
galaxy lensing by galaxies in the Rod-Sequence Cluster Sur- 
vey, Hoekstra et al. (2004) modeled the projected shapes of 
the haloes as ehaio = A eiight , where eught is the ellipticity of 
the image of the luminous galaxy within the halo. Here A = 1 
indicates that the projected shapes of the haloes are identi- 
cal to the shapes of the galaxy images, and A = indicates 
that the haloes axe perfectly circulax in projection on the 
sky. From their analysis, Hoekstra et al. (2004) concluded 
that the haloes of their galaxies were somewhat rounder than 
the images of the galaxies: A = O.TTIq'ji- Using the same 
parametrization of the relationship between the ellipticities 
of the haloes and the images of the galaxies, Mandelbaum et 
al. (2006b) found A = 0.1 ± 0.06 for red SDSS lens galaxies 
and A = — 0.8±0.4 for blue SDSS lens galaxies. Here the neg- 
ative sign indicates an apparent anti-alignment of mass and 
light for blue SDSS lens galaxies. Finally, Parker et al. (2007) 
computed the galaxy-galaxy lensing signal using data from 
the Canada-France-Hawaii Telescope Legacy Survey. When 
Parker et al. (2007) averaged the signal over all lens galax- 
ies, they found a weak (2a) preference for the haloes of the 
lens galaxies to be non-spherical with a projected ellipticity 
of ~ 0.3. When Parker et al. (2007) restricted their anal- 



ysis to elliptical galaxies, the mean halo ellipticity and the 
significance of the result was found to increase somewhat. 

Here we construct a series of Monte Carlo simulations 
in order to explore the theory of weak galaxy-galaxy lens- 
ing by non-spherical dark matter haloes. Using those sim- 
ulations we demonstrate that, in practice, it is challenging 
to interpret the results of an observational effort to detect 
anisotropic galaxy-galaxy lensing. This is because, in gen- 
eral, the observed signature of anisotropic galaxy-galaxy 
lensing is strongly affected by the fact that the central, 
"lens" galaxies have, themselves been weakly lensed. As a re- 
sult, the observed symmetry axes of the central, lens galaxies 
differ from their intrinsic symmetry axes. In our work be- 
low we pay particular attention to the effects of multiple 
weak deflections on the galaxy-galaxy lensing signal. As was 
first pointed out by BBS, galaxy-galaxy lensing is inherently 
a multiple deflection problem. That is, it is common for a 
source galaxy located at redshift Zs to be weakly lensed by 
a galaxy located at zn < Zs- Oftentimes these two galaxies 
are then subsequently lensed by another galaxy at redshift 
zi2 < zii- In other words, the galaxy at zn serves simultar- 
neously as a lens for the galaxy at Zs and a source for the 
galaxy at zi2. In addition, the galaxy at Zs is lensed by two 
different foreground galaxies. Neglecting such multiple de- 
flections when modeling an observed galaxy-galaxy lensing 
signal will give rise to incorrect conclusions about the un- 
derlying properties of the haloes of the lens galaxies. For a 
detailed discussion of the frequency and relative strengths 
of multiple deflections in a deep galaxy-galaxy lensing data 
set, the reader is referred to Brainerd (2010). 

Below, the haloes of the lens galaxies will be modeled as 
truncated singular isothermal ellipsoids. This choice is mo- 
tivated by two considerations. Firstly, the singular isother- 
mal ellipsoid gives rise to a gravitational lensing sheax that 
can be computed analytically (e.g., Kormann, Schneider & 
Bartelmann 1994). Secondly, at the present time the obser- 
vational galaxy-galaxy lensing data are not of sufficiently 
high quality to allow one to distinguish between singular 
isothermal ellipsoid haloes and those that are triaxial CDM 
haloes. 

The outline of the paper is as follows. In Section 2 we 
present the basic theory of gravitational lensing by singular 
isothermal ellipsoids and we introduce a shorthand nota^ 
tion that we will use throughout the paper. In Section 3 
we outline the construction of Monte Carlo simulations of 
galaxy-galaxy lensing by non-spherical haloes, where the lo- 
cations and apparent magnitudes of the Monte Carlo galax- 
ies are taken from a large observational data set. In Sec- 
tion 4 we present the signature of galaxy-galaxy lensing by 
non-spherical haloes that one should expect to obtain from 
a realistic observational data set. In Section 5 we explore 
the effects of galaxy-galaxy lensing on the images of relar 
tively nearby galaxies (i.e., galaxies that are ordinarily be 
considered to bo "lenses" but arc not always considered to 
be "sources"). In Section 6 we construct a second suite of 
Monte Carlo simulations in order to determine the effect 
of uniltiplo weak deflections on observations of anisotropic 
galaxy-galaxy lensing. In Section 7 we demonstrate that the 
effects of lensing of foreground galaxies on the observed sig- 
nature of anisotropic galaxy-galaxy lensing cannot be elimi- 
nated by selective rejection of either lens or source galaxies. 
We summarize our results and present our conclusions in 
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Section 8. Throughout, we will refer to the weak lensing of a 
background galaxy by a single foreground galaxy as a "de- 
flection", and we will adopt a flat A-dominated cosmology 
with Ho = 70 km sec"^ Mpc"\ Qmo = 0.3 and f^Ao = 0.7. 



2 TRUNCATED SINGULAR ISOTHERMAL 
ELLIPSOID LENSES 

Let us assume that the dark matter haloes of large, lumi- 
nous galaxies may be fairly represented as truncated singular 
isothermal ellipsoids. Since we are only concerned with the 
weak lensing regime, the adoption of a halo model that is 
singular (as opposed to a model with a finite density core) 
will have no effect on our results below. Following Kormann 
et al. (1994), the surface mass densities of the dark matter 
haloes are given by 



E(p) 



2G 



1 



(1) 



where a„ is the line of sight velocity dispersion, / is the 

axis ratio of the mass distribution as projected on the sky 
(0 < / ^ 1), xt is the truncation radius, G is Newton's con- 
stant, and p is a generalized elliptical radius defined such 
that = xl + px'i- Here xi and X2 are Cartesian coor- 
dinates measured, respectively, along the minor and major 
axes of the projected mass distribution of the halo. In the 
limiting case of a round lens (i.e., / — ^ 1), the total mass of 
the halo becomes 



Mtot = 



(2) 



The convergence (k) and shear (7 = 71 + 172) are the 
characteristic properties of a gravitational lens. In the case 
of truncated singular isothermal ellipsoid lenses, the conver- 
gence is given by 
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is the critical surface mass den- 



where Ec = j-,^ ^ 

sity, Di is the angular diameter distance of the lens, is the 
angular diameter distance of the source and Du is the angu- 
lar diameter distance between the lens and the source. The 

real and imaginary components of the shear can be obtained 
straightforwardly from equations (63abc) of Kormann ct al. 
(1994): 
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where 




^ _ ■n - - (i - /-)(/'- + ■'7) + ^/-'w + -n 

~ U^r^ - 2/2(1 - P)x^t{^l - xl) + (1 - Pfxt] 

(e.g., Wright 2002). Again, x\ and X2 are Cartesian coordi- 
nates measured along the minor and major axes of the lens, 
respectively. In order to maintain consistency with the no- 
tation of Kormann et al. (1994), here we have used a polar 
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Figure 1. Illustration of our notation, 7+ and 7 (see text). 

coordinate s ystem, c entred on the lens, with radial coordi- 
nate r = 1^ xl + X2 and polar angle, ip, defined such that 
Xl = r cos if and X2 = r sin ip. 

It is clear from equations (4) through (6) that, unlike 
a circularly symmetric lens for which the magnitude of the 
shear depends upon the angular distance from the lens cen- 
tre but not the azimuthal coordinate of the source, the shear 
due to an elliptical lens is a function of both the angular 
distance from the lens centre as well as the azimuthal co- 
ordinate of the source. At a given angular distance, 9, from 
the centre of an elliptical lens, the magnitude of the shear 
is greatest for sources located nearest the major axis of the 
lens and least for sources located nearest the minor axis of 
the lens. Hence, within a given radial anrmlus that is centred 
on the elliptical lens, sources whose azimuthal coordinates, 
(p, place them within ±45° of the major axis of the lens 
will experience a greater mean shear than sources whose az- 
imuthal coordinates, ip, place them within ±45° of the minor 
axis of the lens. As a shorthand notation, we will refer to the 
magnitude of the mean shear experienced by sources whose 
azimuthal coordinates place them within ±45° of the minor 
axis of an elliptical lens as 7" . Similarly, we will refer to the 
magnitude of the mean shear experienced by sources whose 
azimuthal coordinates place them within ±45° of the major 
axis of an elliptical lens as 7+ (see Figure 1). 



3 MONTE CARLO SIMULATIONS OF 
GALAXY-GALAXY LENSING IN THE 
BTC40 SURVEY 

To quantify the effects of non-spherical dark matter haloes 

on the galaxy-galaxy lensing signal, we construct a scries of 
Monte Carlo simulations. As a starting point for our sim- 
ulations we use a set of modestly-deep, wide-field 7-band 
images that were were generously donated to us by Emilio 
Falco. The images were obtained as part of the BTC40 sur- 
vey (Monier et al. 2002), which was carried out using the Big 
Throughput Camera (ETC, Tyson et al. 1992; Wittman et 
al. 1998) on the 4m Blanco telescope at the Cerro-Tololo In- 
teramerican Observatory. The 7-band exposures consist of 
150 seconds per pointing, and each individual image covers 
an area of order one-quarter of a square degree. We selected 
a total of 13.8 deg^ of imaging data from the survey for our 
work, rejecting images that were obtained during poor pho- 
tometric conditions or which exhibited poor tracking or poor 
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focus. The data were calibrated, flat-ficldcd, and de-fringed 
as described in Monier et al. (2002). Object catalogs were 
created from the reduced /-band data using the SExtractor 
package (Bcrtin & Arnouts 1996). 

Additional details regarding the quality of the imaging, 
star-galax;y separation, cosmic ray rejection, point-spread 
function correction, and masking of cosmetic defects (e.g., 
large stellar blooms, diffraction spikes) will be presented in 
a companion paper (Howell & Brainerd, in preparation). In 
the companion paper wc will also present an analysis of the 
observed gala:xy-galaxy lensing signal in this data set. For 
the purposes of our present study, we are simply interested 
in using the BTC40 galaxies as the framework for a set of 
Monte Carlo simulations of galaxy-galaxy lensing by non- 
spherical haloes. That is, here wc will address the follow- 
ing question: Given a data set like that obtained from the 
BTC40, what should one expect to observe for the galaxy- 
galaxy lensing signal if the dark matter haloes of the galaxies 
are non-spherical? The information from the BTC40 images 
that we use here consists solely of the centroids of the galax- 
ies and their 7-band apparent magnitudes. These, along with 
other quantities, are used as input parameters for our Monte 
Carlo simulations. Also, in order to ultimately match the 
data that will be presented in our companion paper, here 
we use only BTC40 galaxies with 18 Iab ^ 22.5. While 
the completeness limit of the data is somewhat fainter than 
Iab = 22.5, in practice the BTC40 galaxies with Iab > 22.5 
are too small for accurate shape determinations. 

The observed shapes of the BTC40 galaxies have been 
affected by the presence of a spatially-varying anisotropic 
point spread function. Because of this, and because of the 
fact that shape determinations become increasingly noisy at 
faint flux levels, wc do not use the observed shapes of the 
BTC40 galaxies in our Monte Carlo simulations. Instead, 
in order to describe the shape of the luminous galaxy, each 
Monte Carlo galaxy is assigned an intrinsic image ellipticity, 
ein = {a — b)/{a + b), that is drawn from the probability 
distribution derived by Ebbels (1998) from 94 archival HST 
field survey images: 



P{t) = ^rexp[-(r/0.036)"-''^]. 



(7) 



Here r = (a^ — 6^)/(2a6), ^ is a normalising constant, and a 
and b are, respectively, the semi-major and semi-minor axes 
of the intrinsic image ellipses. 

We assume that the projected shapes of the haloes of 
the BTC40 galaxies are elliptical but, unlike Hoekstra et al. 
(2004) and Mandelbaum et al. (2006b), we do not assume 
that there is a linear relationship between the shape of the 
luminous galaxy and the shape of its projected dark matter 
halo. While the assumption eimio = Aeught may have some 
validity for elliptical galaxies, it is definitely false for disk 
galaxies (which make up a substantial fraction of the lens 
population). Agustsson & Brainerd (2006) showed that the 
observed ellipticitics of disk galaxies embedded within CDM 
haloes were largely uncorrelated with the ellipticitics of their 
projected haloes (see their Figure 6). This is due to the fact 
that one always views a random projection of the dark mat- 
ter halo on the sky. Therefore, a high inclination angle for the 
disk (which maximises the ellipticity of the luminous galaxy 
image) does not, in general, correlate with a projection that 
maximises the projected ellipticity of the halo. In order to 
assign projected axis ratios, /, to the haloes of our Monte 
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projected halo axis ratio, f 

Figure 2. Distribution of projected axis ratios, /, for the Monte 
Carlo dark matter haloes (following the results of Agustsson & 
Brainerd 2006). 



Carlo galaxies, then, wc use the probability distribution ob- 
tained by Agustsson & Brainerd (2006) for the projected 
axis ratios of CDM galaxy haloes. The halo of each galaxy 
in our simulations is therefore assigned a value of / that is 
drawn at random from this distribution (see Figure 2). 

Next wc must make a choice as to how to orient the 
luminous galaxies within their dark matter haloes. The only 
symmetry axes that can be used in an observational data set 
to detect anisotropic galaxy-galaxy lensing are, of course, the 
symmetry axes of the luminous galaxies themselves. If mass 
and light arc not reasonably well aligned within the lens 
galaxies, a detection of anisotropic galaxy-galaxy lensing is 
hopeless since we cannot directly observe the orientations of 
the symmetry axes of the dark matter haloes. Therefore, in 
our simulations we will assume that the intrinsic symmetry 
axes of the luminous galaxies and their dark matter haloes 
are aligned with each other. This assumption maximises the 
degree of anisotropy in the galaxy-galaxy lensing signal that 
one should expect to see and it presents a best case scenario 
for detecting the effect. 

Neither spectroscopic redshifts nor photometric red- 
shifts axe available for the BTC40 galaxies. Therefore, we 
must assign redshifts to the galaxies in order to carry out 
our Monte Carlo simulations. Following the prescriptions of 
BBS and Wright (2002), we adopt a redshift distribution of 
the form 



P{z\Iab) 



iSz^e^p [-{z/zof] 
T(3/I3)zl 



(8) 



Taking /3 = 1.5 yields good agreement with the redshift 
surveys of LeFevre et al. (1996) and LeFevre et al. (2004), 
and we then have 



zo 



0.8 [0.86 + 0.15{Iab - 23.35)] , 



(9) 



Lastly, we must assign velocity dispersions and trunca- 
tion radii to the haloes of each of the Monte Carlo galax- 
ies. To do this, wc assume that the galaxies follow a Fabcr- 
Jackson or TuUy-Fisher type of relationship and have con- 
stant mass-to-light ratio (see, e.g., BBS). The velocity dis- 
persion, cr„, of the halo of a galaxy with luminosity, L, is 
then given by 



at ~ VLV 



1/4 



(10) 
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where is the velocity dispersion of the halo of an L* 
galaxy. The truncation radius, xt, of the halo of a galaxy 
with luminosity, L, is given by 




where Xt is the truncation radius of the halo of an L* galaxy. 
The luminosity of each Monte Carlo galaxy is obtained from 
its observed J-band apparent magnitude and the redshift, z, 
that was assigned to the galaxy based on equation ((8| above. 
Accounting for the K-correction, we have 

where a = ——^^7—^ (e-g-i BBS). For simplicity, we take 
a = 0.42, which is the mean slope of the spectral energy dis- 
tribution between the Johnson _R-band and B-band from the 
Caltech Faint Galaxy Redshift Survey (Cohen et al. 1999ab). 
For each Monte Carlo simulation then: 

• a pair of characteristic parameters, {a*,Xt), are 
adopted for the haloes of L* galaxies 

• each luminous galaxy is assigned its observed location 
on the image, as well as its observed /-band apparent mag- 
nitude 

• the image of each luminous galaxy is assigned an intrin- 
sic shape, eiii = (a — b)/{a + b), using equation ^ and its 
dark matter halo is assigned an axis ratio, /, drawn from 
the projected halo shapes in Agustsson & Brainerd (2006) 

• each luminous galax;y is assigned a random intrinsic po- 
sition angle, (pin, (i.e., we assume that in the absence of grav- 
itational lensing the galaxy images are uncorrelated) and, 
since we also assume that mass and light are aligned in pro- 
jection on the sky, the projected halo is assigned a position 
angle identical to the position angle of the unlensed lumi- 
nous galaxy 

• each galaxy is assigned a redshift, z, using equation ((HJ, 
and its luminosity relative to L* is obtained using equation 

m 

• the dark matter halo of each galaxy is assigned a ve- 
locity dispersion, cr„, and truncation radius, xt, based upon 
the luminosity of the galaxy within the halo and the scaling 
relations of equations (I10|l and (1111) 

Each Monte Carlo simulation then proceeds by computing 
the weak lensing shear, 7, that is induced as light rays ema- 
nating from distant galaxies encounter the gravitational po- 
tentials of foreground galaxies. As we will see below, most 
of the distant galaxies with redshift Zi are lensed by numer- 
ous foreground galaxies with redshifts Zj < Zi. We define 
the intrinsic (unlensed) shape of each luminous Monte Carlo 
galaxy to be 

Xin = tin e^***'" (13) 

where ein in the intrinsic (unlensed) ellipticity of the galaxy 
image and <f)in is the intrinsic (unlensed) position angle. 
Since we are dealing with the weak lensing regime, all lensing 
events may be considered to be independent (e.g., Bartel- 
mann & Schneider 2001) and the final image shape of each 
lensed galaxy is given by 

Xf = Xin + S^^-" 7j = Xin + Xnot (14) 

where yj is the shear induced by foreground lens galaxy, j, 
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Figure 3. Redshift distributions adopted for the "bright" fore- 
ground (18 ?J Iab ^ 20) and "faint" background (20 < Iab 'Si 
22.5) galaxies in the BTC40 sample. 



and Xnct is the net shear due to all foreground lenses. The 
real (71) and imaginary (72) components of the shear are 
given by equations (4) through (6) above. 

Computation of the net shear for each of the galaxies 
due to literally all potential foreground lens galaxies is ex- 
tremely time-consuming and, from a practical standpoint, 
is unnecessary since foreground lenses that induce negligi- 
ble shear (say, ~ 10~^) can be neglected in compari- 
son to foreground lenses that induce substantial shear (say, 
7j > 0.005). From Brainerd (2010), we know that source 
galaxies with a median redshift Zs = 0.96 that have been 
lensed by a population of foreground galaxies with Zi = 0.55 
experience little shear due to lenses that are located at pro- 
jected radii 6 > 60". Scaling to the BTC40 gala:xies, we 
find that for 6 > 100" the contribution to the net galaxy- 
galaxy lensing shear will be negligible. Hence, in our Monte 
Carlo simulations we compute the net shear experienced by 
each BTC40 galaxy due to all foreground galaxies that are 
located within a projected radius of 100", and we do not in- 
clude any contribution to the net shear from lenses located 
at projected radii > 100". 

In the next section we will analyse the output of our 
Monte Carlo simulations in a manner that is similar to the 
way in which an observational galaxy-galaxy lensing data set 
is analysed. In the case that neither spectroscopic nor pho- 
tometric redshifts are available for an observational data set 
(as is the case for the BTC40 data), one can make only a 
crude distinction between "foreground" and "background" 
galaxies using apparent magnitudes. From the probability 
distribution above, we know that, on average, galaxies with 
bright apparent magnitudes tend to be located at lower red- 
shifts than galaxies with faint apparent magnitudes (though 
there is certainly a good deal of overlap). If we consider 
galaxies with 18 ^ Iab ^ 20, we find a median redshift of 
Zmcd ~ 0.29 for our BTC40 sample. If we consider galax- 
ies with 20 < Iab ^ 22.5, we find a mean redshift of 
2mcd ~ 0.61 for our BTC40 sample. These magnitude cuts 
therefore yield a rough division of our BTC40 sample into 
38,879 "bright" foreground (18 sC Iab ^ 20) objects and 
225,518 "faint" background (20 < Iab < 22.5) objects. The 
redshift distributions adopted for these objects are shown 
Figure 3. 

It is important to remember that in the Monte Carlo 
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Figure 4. Results of one Monte Carlo simulation for one CCD 
frame from the BTC40 data. Brown dots show the locations of 
bright, foreground galaxies (18 ^ Iab ^ 20) and blue dots show 
the locations of faint, background galaxies (20 < Iab ^ 22.5). 
Lines indicate the direction on the sky of objects that have lensed 
a given galaxy. Blue lines indicate that the lens is a faint galaxy; 
brown lines indicate that the lens is a bright galaxy. Note that 
virtually all bright, foreground galaxies have been lensed multiple 
times. 

simulations, all galaxies with redshifts Zi have been lensed 
by all other galaxies with redshifts Zj < Zi that are located 
within a projected radius 6 — 100" of the galaxy at redshift 
Zi. We illustrate this in Figure 4, where we show an exam- 
ple of a Monte Carlo simulation. The image corresponds to 
a single CCD frame from the BTC40 survey (0.25 deg x 
0.25 deg), where the dots indicate the locations of galax- 
ies with magnitudes 18 ^ Iab ^ 22.5 and the lines indi- 
cate the directions on the sky of objects that have lensed 
the galaxies. That is, each line indicates the presence of 
a lens-source pair in the data; however, for clarity of the 
figure we do not extend the lines to connect every source 
galaxy directly to all of its lenses. Brown dots show the lo- 
cations of bright, foreground galaxies (18 55 Iab ^ 20) and 
blue dots show the locations of faint, background galaxies 
(20 < Iab ^ 22.5). Blue lines indicate that the lens is a faint, 
background galaxy. Brown lines indicate that the lens is a 
bright, foreground galaxy. Therefore, a blue line originating 
from a blue dot indicates that a faint, background galaxy has 
been lensed by another faint, background galaxy. Similarly, 
a brown line originating from a blue dot indicates that a 
faint, background galaxy has been lensed by a bright, fore- 
ground galaxy. Importantly, Figure 4 shows that virtually 
all of the bright, foreground galaxies (the brown dots) have, 
themselves, been lensed multiple times. Most of the bright, 
foreground galaxies have been lensed by other bright galax- 
ies but they are occasionally lensed by a faint galaxy (due to 
the overlapping redshift distributions of these objects). The 
majority of the lenses turn out to be faint galaxies simply 
because there are ~ 6 times as many faint galaxies per unit 
area as there are bright galaxies. 

It is clear from Figure 4 that the vast majority of the 



galaxies have been lensed by more than one galaxy; that is, 
multiple deflections are common for all of the galaxies in the 
Monte Carlo simulations. This statement is true indepen- 
dent of the values of the characteristic parameters adopted 
for the haloes of L* galaxies, {a'^,Xt); however, the relative 
strengths of the individual deflections and their net effect 
on Xf for each galaxy will, of course, be a strong function of 
the values of the characteristic parameters that are adopted 
(see, e.g., Brainerd 2010). 



4 SIGNATURE OF ANISOTROPIC 
GALAXY-GALAXY LENSING 

Here we use the output of the BTC40 Monte Carlo simu- 
lations to compute the dependence of the mean tangential 
shear of the faint (20 < Iab ^ 22.5) galaxies using the 
bright (18 55 Iab ^ 20) galaxies as the centres for the calcu- 
lation. That is, we compute the signature of galaxy-galaxy 
lensing in the same way as is done for an observational data 
set in which apparent magnitude is used as the sole dis- 
criminator between "foreground" and "background" objects. 
We separately compute j'^iO) and 7~(^), and we show the 
results in Figures 5-7 for three different sets of character- 
istic parameters that were chosen to represent the haloes 
of L* galaxies. Since galaxy-galaxy lensing is relatively in- 
sensitive to the radii of the dark matter haloes (e.g., BBS; 
Hoekstra et al. 2004; Brainerd 2010), we adopt a value of 
Xt = 100 kpc for the haloes of L* galaxies. In Figures 5-7 
we then vary the characteristic velocity dispersion, adopting 
values of = 100 km sec~^ (Figure 5), = 150 km sec~^ 
(Figure 6), and cr* — 200 km sec~^ (Figure 7). 

In order compute 7^ (9) and 7" (9) we must flrst deflne 
what we mean by the symmetry axes of the bright galax- 
ies. Implicit in our definition of 7"*" and 7" is that the in- 
trinsic (unlensed) symmetry axes are the symmetry axes of 
the projected dark matter halo (e.g.. Figure 1). However, 
observers are not blessed with "dark matter glasses" that 
allow us to see the intrinsic symmetry axes. Rather, in an 
observational data set, we must take the symmetry axes of 
the bright galaxies to be their observed symmetry axes, not 
their intrinsic symmetry axes. This is an important distinc- 
tion since the observed symmetry axes of the bright galaxies 
may have been altered due to weak lensing by foreground 
galaxies; see, e.g.. Figure 4. 

In Figures 5-7 we compute 7"^ (9) and 7^ (9) using both 
the observed symmetry axes (circles) and the intrinsic sym- 
metry axes (crosses) of the bright, central galaxies. That is, 
the circles indicate the functions that we would expect to 
measure in an observational data set, while the crosses in- 
dicate the functions that we would obtain if we were able 
to observe the intrinsic (unlensed) symmetry axes of the 
bright, central galaxies. In the case of very low mass lenses 
(Figure 5), there is relatively little difference between the 
tangential shears that result from using the observed sym- 
metry axes of the bright centres and those that result from 
using the intrinsic symmetry axes. However, for moderate to 
high mass lenses (Figures 6 and 7), it is clear that over most 
scales there is a systematic difference between the two calcu- 
lations. In particular, over most scales the observed values of 
"/^{d) in Figures 6 and 7 are systematically lower than than 
the values that are obtained by using the intrinsic symmetry 
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Figure 5. Observed functions, -f'^{6) and 7^(6), from the BTC40 
Monte Carlo simulations. Results show the mean of 15 indepen- 
dent realisations of the complete BTC40 data set. All source 
galaxies that are located within ±45° of the symmetry Eixes of 
the bright, foreground centres are used in the calculations. Dif- 
ferent point types indicate different definitions of the symmetry 
axes of the bright, foreground centres (circles: observed symmetry 
axes after lensing, crosses; intrinsic symmetry axes). Error bars 
arc omitted because they are smaller than the data points. Here 
the characteristic parameters for the haloes of L* galaxies are 
a* = 100 km sec"! and x* = 100 h'^ kpc. 



axes of the bright centres. Conversely, over most scales the 
observed values of 7" (9) in Figures 6 and 7 are systemati- 
cally higher than than the values that are obtained by using 
the intrinsic symmetry axes of the bright centres. 

Shown in Figure 8 is the ratio of the mean tangential 
shears, "f'^ (0) / 'y' (9) , for our three halo models. From this 
figure, we see that for low mass haloes {al = 100 km sec~^) 
the ratio of the mean tangential shears is slightly lower on 
average when the shear is measured with respect to the ob- 
served symmetry axes of the bright centres than when it 
is measured with respect to the intrinsic symmetry axes. 
However, to within the error bars, the two functions for- 
mally agree. For moderate (cr* = 150 km sec~^) to high 
(ct* = 200 km sec^"*^) mass lenses there is a substantial 
suppression of 7^(^)/7~(^) when the observed symmetry 
axes of the bright centres are used compared to what one 
would obtain using the intrinsic symmetry axes. In the case 
of cr* = 150 km sec~^ there is little to no anisotropy ap- 
parent on scales 9 > 20". That is, the observed signature 
of anisotropic galaxy-galaxy lensing by haloes of moderate 
mass is largely consistent with isotropic galaxy-galaxy lens- 
ing on scales 9 > 20". In the case of cr^ = 200 km sec~^, 
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Figure 8. Ratio of the mean tangential shears, 7^(f)/7" (9), for 
three halo models adopted for the BTC40 galaxies. In all eases 
the truncation radius of the haloes of L* galaxies is taken to be 
X* = 100 h^^ kpe. Here all source galaxies that are located within 
±45° of the symmetry axes of the bright, foreground centres are 
used in the calculation. Different point types indicate different 
definitions of the symmetry axes of the bright, foreground centres 
(circles: observed symmetry axes after lensing, diamonds: intrinsic 
symmetry axes). Error bars are omitted when they are compara- 
ble to or smaller than the data points. Top: a* = 100 km sec~^. 
Middle: at = 150 km sec~^. Bottom: a*. = 200 km sec~^. 



the observed function is actually reversed from the expected 
function (i.e., ^+{6) < 7^(6^)) on scales 20" < 6 < 70". 
That is, the observed signature of anisotropic gala^xy-galEixy 
lensing by high mass haloes could lead one to think (falsely) 
that mass and light are anti-aligned within the galaxies. 

Figure 8 demonstrates, then, that an observation of 
7"'"(0) = 7~(5) is not categoric proof that the haloes of 
the lenses are spherically- symmetric since the haloes of our 
Monte Carlo galaxies are non-spherical. In addition, an ob- 
servation of 7^(^) < 1~ {&) is not categoric proof that mass 
and light are anti-aligned within the lens galaxies since the 
intrinsic symmetry axes of the luminous galaxies in the 
Monte Carlo simulations were taken to be aligned with the 
symmetry axes of their projected dark matter haloes. Since 
we have allowed the bright, central galaxies that we have 
used to compute 7^ {9) and 7" {9) to be lensed by foreground 
galaxies, the circles in Figure 8 show the actual signature of 
galaxy-galaxy lensing by non-spherical dark matter haloes 
that one should expect to see in an observational data set 
(i.e., a data set in which the galaxies are all broadly dis- 
tributed in redshift space). Contrary to the usual expecta- 
tion that 7"'"(^) should exceed 7^(^) over a wide range of 
angular scales, Figure 8 shows that this is unlikely to be the 



case unless the haloes of L* galaxies have particularly low 
characteristic velocity dispersions {a^ — 100 km sec^^). In 
Section 6 we will demonstrate that the results shown in Fig- 
ure 8 are caused primarily by the fact the observed symme- 
try axes of the bright, foreground centres have been altered 
from their intrinsic symmetry axes by weak lensing. 

Figures 9 and 10 show schematic illustrations of what 
can occur in the situation that a given lens-source pair is, 
itself, lensed by a foreground mass. Consider an intrinsically 
circular source galaxy that is located near the major axis 
of an elliptical lens galaxy. That is, the lens-source pair is 
in what one might call the "7^ configuration". After being 
sheared by the elliptical lens, the image of the source galaxy 
is an ellipse with the major axis of its image oriented tan- 
gentially with respect to the major axis of the elliptical lens. 
This is illustrated by the blue ellipse in the top and bottom 
panels of Figure 9. Now consider the effect on the image of 
the intrinsically circular source if a large mass (i.e., another 
galaxy) is placed in the foreground of the original lens-source 
pair. If the additional foreground mass is located along the 
line that connects the centroids of the original lens-source 
pair, the net result for the image of the intrinsically-circular 
source is that it becomes more elliptical than if it had been 
lensed only once. This is illustrated by the red ellipse in the 
top panel of Figure 9. Now consider the effect if the fore- 
ground mass is placed such that its location is tangential 
to the line that connects the centroids of the original lens- 
source pair. The net result for the image of the intrinsically- 
circular source is that it will be less elliptical than if it had 
been lensed only once. This is illustrated by the red ellipse 
in the bottom panel of Figure 9. We see, therefore, that the 
inclusion of a second lens may either increase or decrease 
the net ellipticity of our distant, circular source over what 
we would have naively expected in the single-deflection case. 

Next let us consider the fact that not only will the in- 
troduction of an additional foreground mass alter the image 
of the distant, circular source galaxy, it will also alter the 
image of the original elliptical lens galaxy. Again, consider 
the original lens-source pair to be in the 7"*" configuration. 
If the foreground mass is placed along the line that connects 
the centroids of the original lens-source pair, the image of 
the elliptical lens becomes rounder than its intrinsic shape 
(i.e., since it is distorted tangentially with respect to the 
location of the foreground mass). In some finite number of 
cases where the intrinsic ellipticity of the image of the ellip- 
tical lens is very small and the shear due to the foreground 
lens is large, the post-lensing image of the elliptical lens may 
even have its observed symmetry axes reversed from its in- 
trinsic symmetry axes. As a result, the distant, intrinsically 
circular source galaxy would appear to be iu the 7^ con- 
figuration when, in fact, it is in the 7"*" configuration (e.g., 
top panel of Figure 10). Should the situation illustrated in 
the top panel of Figure 10 occur, the mean tangential shear 
that one would obtain for sources located close to the 06- 
served minor axes of the elliptical lens will be greater than it 
ought to be. That is, the observed value of 7" is boosted by 
the incorrect inclusion of sources that would have properly 
gone into the calculation of 7+ if one had known the orien- 
tation of the intrinsic symmetry axes of the elliptical lens. 
Of course, this also results in the observed value of 7"*" be- 
ing reduced compared to its true value because the observed 
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Figure 9. Schematic illustration of the net effect on the image of 
a distant, intrinsically-circular source galaxy due to an elliptical 
lens with redshift z;i < Zs and an additional, foreground lens with 
redshift Z12 < zn. The black ellipse indicates the true shape and 
orientation of the elliptical lens. In both panels the elliptical lens 
and the source are intrinsically in the "7+" configuration. Blue 
ellipses: shape of the source after lensing solely by the elliptical 
lens. Red ellipses: shape of the source after being lensed by the 
elliptical lens at zn and the additional foreground lens at 2;2. 
When the foreground lens is located along the direction vector 
that connects the centroids of the elliptical lens and the circular 
source, the final ellipticity of the image of the source is increased 
compared to what it would have been if the source had been lensed 
solely by the elliptical lens. When the foreground lens is located 
tangential to the direction vector that connects the centroids of 
the elliptical lens and the circular source, the final ellipticity of 
the image of the source is reduced compared to what it would 
have been if the source had been lensed solely by the elliptical 
lens. 




Lens 2 (z,2 < z,^) Lens 1 (z,, < z^) Source (z^) 




Figure 10. Schematic illustration of the net effect on (i) the 
image of a distant, intrinsically-circular source galaxy due to an 
elliptical lens with redshift Zn < z^ and an additional, foreground 
lens with redshift Z12 < zn, and (ii) the image of the elliptical 
lens after being lensed by the foreground lens at 2; 2- The black 
ellipse indicates the true shape and orientation of the elliptical 
lens in both panels. Top: elliptical lens and source are intrinsically 
in the "'y+" configuration. Bottom: elliptical lens and source are 
intrinsically in the "7"" configuration. Blue ellipses: shape of the 
image of the source after lensing solely by the elliptical lens. Red 
ellipses: shape of the image of the elliptical lens after being lensed 
by the foreground lens at Z12, as well as the shape of the image of 
the source after being lensed by the elliptical lens at Zn and the 
additional foreground lens at Zi2. This figure illustrates that in the 
limit of small intrinsic ellipticities, the observed symmetry axes 
of the elliptical lens may be reversed from its intrinsic symmetry 
axes. This results in a misclassification of the configuration (7+ or 
7~) of the elliptical lens-source pair, which incorrectly enhances 
measurements of ■y~ and incorrectly suppresses measurements of 



lens-source configuration has been "misclassified" compared 
to its intrinsic configuration. 

Next consider placing the foreground mass along a line 
that is tangential to the line that connects the centroids 
of the original lens-source pair. If the original lens-source 
pair is in the 7^ configuration, the image of the elliptical 
lens will have an increased ellipticity after being lensed by 
the foreground mass and it will not undergo a reversal of 
its symmetry axes (i.e., the lens-source pair remains in the 
7"*" configuration), but the final image of the source will be 
less elliptical than it would have been in the absence of the 
foreground lens. In the case that the lens-source pair is in- 
trinsically in the 7" configuration, however, a finite num- 
ber of lenses with small intrinsic ellipticities may have their 
symmetry axes reversed by lensing due to the foreground 
mass (e.g., bottom panel of Figure 10). Hence, the original 



lens-source pair would appear to be in the 7 configuration, 
when, in fact, it is actually in the y~ configuration. In this 
case, the intrinsically circular source is incorrectly put into 
the calculation of 7^, and its net shear will be rather small 
because: (i) it is, in reality, located near the minor axis of 
the elliptical lens and (ii) after being lensed by the addi- 
tional foreground mass, its image will be rounder than if it 
had been lensed solely by the original elliptical lens. Both of 
these conspire to reduce the observed value of 7^ compared 
to its true value. 

Based upon relatively simple reasoning from Figures 9 
and 10, we may therefore expect that multiple deflections 
(i.e., the presence of more than one foreground lens) could 
lead to a suppression of {d)/j~ (8) compared to what 
would would obtain if one knew the intrinsic symmetry axes 
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Figure 11. Transformation of an ellipse by a scalar shear. Blue: 
original ellipse. Red: transformed ellipse. Note that not only are 
the ellipticity and position angle of the ellipse altered, the original 
major and minor axes of the ellipse are no longer orthogonal. 

of the bright centres. One might hope that in a sufficiently 
large data set the effects of foreground lenses on the original 
elliptical lens-source pair would cancel each other out. How- 
ever, this is not necessarily going to be the case. Galaxies 
span a broad range of redshifts; hence, at fixed angular sep- 
aration from a source, 6, two foreground lens galaxies with 
identical gravitational potentials will have different lensing 
strengths because they are located at different physical dis- 
tances from the source. In the following sections we will ex- 
plore the effect of galaxy-galaxy lensing on the images of 
the bright, foreground objects that are used as centres for 
the computation of the mean tangential shear of the faint, 
background objects. In addition, we will explore the effect of 
weak lensing of the central galaxies on the measured values 
of 7+(6l) and 7^(61). 
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Figure 12. Transformation of an ellipse of intrinsic ellipticity, 
tin, due to a scalar shear of magnitude 7 = 0.01. Prior to being 
sheared, the position angle of the ellipse is ^. The quantity is 
the angle between the major axis of the unsheared source ellipse 
and the direction vector that connects the centroids of the lens 
and source. Top: change in position angle of the ellipse. Bottom: 
ratio of intrinsic and transformed ellipticities. 



5 EFFECT OF GALAXY-GALAXY LENSING 
ON FOREGROUND GALAXY IMAGES 

All of our Monte Carlo galaxies with redshifts Zi have been 
lensed by all other Monte Carlo galaxies with redshifts 
Zj < Zi that are found within a radius of 100" of the 
galaxy at Zi. Since the redshift distribution of the galax- 
ies is broad, this means that many of the bright centres 
(i.e., those BTC40 galaxies with 18 ^ Iab ^ 20) corre- 
spond to galaxies that have, themselves, been lensed (e.g., 
Figure 4). At a given angular separation from a weak galaxy 
lens, galaxy-galaxy lensing may be considered to be a scalar 
shear. That is, although a source galaxy has a finite size on 
the sky, the shear due to a foreground weak galaxy lens is 
effectively constant across the image of the source. Shown 
in Figure 11 is an illustration of the transformation of an 
ellipse (i.e., the intrinsic shape of a galaxy) due to a scalar 
shear. The important things to note from this figure are that 
a scalar shear applied to an ellipse results in a change in the 
ellipticity as well as a change in the position angle. In addi- 
tion, the original major and minor axes of the ellipse are no 
longer orthogonal after the transformation. 

In the limit that the intrinsic ellipticity, tin = (a — 
b)/{a + b), of an ellipse is large compared to the applied 
shear, tin » 7/2, the transformation of the ellipse due to 
a scalar shear can be obtained straightforwardly from equa- 
tions (39) and (40) of Surpi & Harari (1999). Let the ellipse 



have intrinsic axis ratio / = b/a and intrinsic position angle 
tp. If we then take to be the angle between the major axis 
of the unlensed source ellipse and the vector that connects 
the centroids of the lens and the source ellipse, the resulting 
change in position angle of the source ellipse is 

AtA- 7(^1^^ sin 2(0-^) (15) 

and the square of the axis ratio of the transformed source 
ellipse is given by 

(f'f _ .f'-27/'cos2(0-7A) 

1 + 27/2 cos 2(<j!>-^) • ' 

Figure 12 shows the resulting change in position angle and 
ellipticity for ellipses with intrinsic ellipticity 0.02 ^ ein ^ 
0.3 due to a 1% scalar shear (7 = 0.01). Unsurprisingly, the 
smaller is ein, the greater is the change in position angle, 
and the maximum change in the position angle occurs for 
{(f) — ip) ~ 45°. In addition, the smaller is ein, the greater is 
the change in ellipticity. The maximum change in ellipticity 
occurs for the two extreme conditions: {(j)—tp) ~ 0°, resulting 
in an ellipse that is rounder than its intrinsic shape, and 
(0 — V) ~ 90° J resulting in an ellipse that is flatter than its 
intrinsic shape. 
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6 EFFECT OF MULTIPLE WEAK 

DEFLECTIONS ON 7+ (61) AND ^-{6) 

In the previous section we demonstrated that, in the Umit 
of small intrinsic ellipticities, the images of the galaxies that 
correspond to the "bright centres" in our galaxy-galaxy lens- 
ing analysis could be significantly affected by weak lensing 
by foreground galaxies under the right circumstances. That 
is, multiple weak deflections (in which a lens-source pair is 
subsequently lensed by one or more foreground galaxies) will 
affect not only the image of the original source, but they will 
also affect the image of the original lens. In this section wc 
investigate the effect of such multiple weak deflections on 
measurements of 7"'"(0) and j~{0)- Here we wish to sepa- 
rate the effects of weak lensing due to non-spherical dark 
matter haloes from the effects of subsequent multiple weak 
deflections due to the haloes of foreground galaxies. To ad- 
dress this, we construct a set of constrained Monte Carlo 
simulations in which a lens-source pair that contains an el- 
liptical lens is, itself, weakly-lensed by additional galaxies 
whose dark matter haloes are spherically-symmetric. 

This second set of Monte Carlo simulations is con- 
structed as follows. An elliptical lens of fixed mass axis ratio, 
/, is placed at the origin of the coordinate system. Here / is 
the axis ratio of the projected dark matter halo. The ellipti- 
cal lens is assigned a random position angle and a redshift, 
Zei, drawn from the redshift distribution adopted for the 
BTC40 objects with apparent magnitudes 18 ^ Iab ^ 20. 
The elliptical lens is assigned an intrinsic, luminous galaxy 
shape drawn from equation (7) above, and the intrinsic po- 
sition angle of the luminous lens galaxy is taken to be the 
position angle of its dark matter halo. A source galaxy is 
placed randomly along the horizontal axis of the coordinate 
system, such that it lies within a distance of ±30" of the 
elliptical lens. The source galaxy is assigned a fixed redshift, 
Zs = 0.6. 

Next, a circle of radius 60", centred on the elliptical lens, 
is populated with additional galaxy lenses whose haloes are 
taken to be singular isothermal spheres (SIS). The number 
density of the SIS lenses is matched to the observed number 
density of galaxies with apparent magnitudes 18 ^ Iab ^ 20 
in the BTC40 data, and they are assigned random locations 
within the field. Each SIS lens is assigned a redshift based 
upon its apparent magnitude, again drawn from our adopted 
probability distribution. 

The elliptical lens is assigned a velocity dispersion, 
and truncation radius, xt. For simplicity, the SIS lenses are 
assigned a velocity dispersion equal to the velocity disper- 
sion that is assigned to the elliptical lens. Having assigned 
positions, redshifts, and gravitational potentials to all of the 
lenses, then, the net shear experienced by the source is com- 
puted as the sum of the individual shears due to all lenses 
(elliptical and SIS) with redshifts zi < 0.6. In addition, the 
final image shape of the lens galaxy residing within the ellip- 
tical dark matter halo is computed using the net shear due 
to all foreground SIS lenses (i.e., SIS lenses with zt < Zei). 

The above procedure is repeated 20 million times for a 
fixed axis ratio, /, of the halo of the elliptical lens, fixed val- 
ues of CT„ and Xt, and fixed source redshift, Zs = 0.6. For each 
new Monte Carlo realisation, a new intrinsic position angle 
and a new intrinsic image shape for the luminous galaxy 
within the elliptical lens halo are generated. That is, with 
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Figure 13. Ratio of mean tangential shears, 7''"('?)/7~ (9), for 
our suite of constrained Monte Carlo simulations. Here the halo 
of the central, elliptical lens galaxy has ehalo = O-l- The func- 
tion is measured relative to the symmetry axes of the central lens 
galaxy. Different panels show results for various values of the ve- 
locity dispersion and truncation radius of the central lens galaxy. 
All source galaxies located within ±45° of the syimnctry axes of 
the central lens galaxy arc included in the calculation. Error bars 
are omitted because they arc comparable to or smaller than the 
data points. Circles; Source galaxies have been Icnscd solely by 
the central, elliptical lens. Symmetry axes used in the calculation 
of l^{Q)/^~{0) are the intrinsic symmetry axes of the central 
lens. Crosses: Source galaxies have been lensed by the central, 
elliptical lens as well as all foreground SIS lenses. Symmetry axes 
used in the calculation of 7~''('')/7~(^^) are the intrinsic symme- 
try axes of the central lens. Squares: "Observed" signal. Source 
galaxies have been lensed by the central, elliptical lens as well 
as all foreground SIS lenses. Image of the central, elliptical lens 
has also been lensed by all foreground SIS lenses. Symmetry axes 
used in the calculation of 'y'^ (9) / ^~ (Q) arc the final, observed 
symmetry axes of the central, elliptical lenses after lensing by the 
foreground SIS lenses. 



each new realisation we randomly "spin" the elliptical lens 
and we assign its luminous galaxy a new intrinsic elliptic- 
ity. In addition, the elliptical lens is assigned a new redshift 
in each new realisation. Thus, after many realisations, the 
redshifts of the elliptical lenses will span the entire range 
of redshift space that was adopted for BTC40 galaxies with 
18 < Iab ^ 20. For each new Monte Carlo realisation a new 
location for the source along the horizontal axis is gener- 
ated, and a new suite of SIS lenses is laid down (including 
new redshifts and new locations within the 60" circle). After 
the source at Zs = 0.6 and the elliptical lens at the origin 
have been lensed by all foreground galaxies, the net shear 
for both the source and the luminous galaxy within the el- 
liptical halo are computed in each individual Monte Caxlo 
realisation. The mean tangential shear for the sources, 7"*" 
and 7~ , is then computed by taking the elliptical lenses as 
the centres for the calculation. 

Shown in Figures 13-15 is the function {0)/j~ (0), ob- 
tained by computing the mean tangential shear for sources 
located within ±45° of the symmetry axes of the ellipti- 
cal lens. Results from a range of halo parameters ((t„ = 
100 km sec"\ 150 km sec"\ 200 km scc"^ xt = 50 kpc, 
100 h-^ kpc, 200 h-^ kpc) are shown in the different panels. 
Figure 13 shows results for central elliptical lenses in which 
the ellipticity of the dark matter halo is ehaio =0.1 (corre- 
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sponding to a projected mass axis ratio / = 0.82). Figure 14 
shows results for central elliptical lenses in which the ellip- 
ticity of the dark matter halo is ehaio ~ 0.3 (corresponding 
to a projected mass axis ratio of f — 0.54). Figure 15 shows 
results for central elliptical lenses in which the cllipticity 
of the dark matter halo is ehaio = 0.5 (corresponding to a 
projected mass axis ratio of / = 0.33). 

The circles in Figures 13-15 show f'^iO)/'y~{d) for the 
case that the source galaxies are lensed solely by the central 
elliptical lens. The symmetry axes used for the calculation 
are the intrinsic symmetry cixes of the central elliptical lens. 
In other words, the circles show the simplest expected re- 
sult: all source galaxies arc Iciiscd by only one foreground 
galaxy, the lens has a non-spherical dark matter halo, and 
the symmetry axes of the lens galaxy are its intrinsic sym- 
metry axes. The crosses in Figures 13-15 show 'f'^ {9) (d) 
for the case that the source galaxies are lensed by both the 
central elliptical lens, as well as all foreground SIS lenses. 
The symmetry axes used for the calculation are the intrin- 
sic symmetry axes of the central elliptical lens. Comparing 
the crosses to the circles we find that the introduction of 
foreground SIS lenses does little to affect the ratio of the 
tangential shears when the intrinsic symmetry axes of the 
central elliptical lens are used for the calculation. 

The squares in Figures 13-15 show 'y^{0)/^~{9) for the 
case that the source galaxies are lensed by both the cen- 
tral elliptical lens, as well as all foreground SIS lenses. In 
addition, the central, elliptical lens has been lensed by all 
foreground SIS lenses. Here the symmetry axes used for the 
calculation are the observed symmetry axes of the central el- 
liptical lens (i.e., the symmetry axes after lensing by the fore- 
ground SIS lenses). From Figures 13-15, the degree to which 
the observed function, 7^(^)/7 (^), is suppressed compared 
to what one would obtain using the intrinsic symmetry axes 
of the elliptical lens is a function of the velocity dispersion 
that is adopted. The lower is the velocity dispersion of the 
lenses, the less the observed function is suppressed. This is 
due to the fact that the frequency and strength of the mul- 
tiple weak deflections are lower for lenses with low veloc- 
ity dispersions than for lenses with high velocity dispersions 
(sec, e.g., Brainerd 2010). In contrast, the cllipticity of the 
projected dark matter halo of the central elliptical lens has 
relatively little effect on the degree to which the observed 
function, ')'^ {0) / 'y~ (0) is suppressed. (Note that the vertical 
scales in Figures 13-15 are very different from each other due 
to the fact that the more elliptical is the central elliptical 
lens, the greater is the anisotropy that it induces.) 



7 JUDICIOUS REJECTION OF LENSES AND 
SOURCES? 

In an attempt to "inoculate" one's analysis of galaxy-galaxy 
lensing by non-spherical haloes against the above effects, one 
might consider simply rejecting bright centres with small 
image ellipticitics from the calculation oi (6) and ^~{6). 
That is, one could hope to avoid the extreme situation where 
a lens-source pair that is truly in the "7+" configuration is 
swapped to the "7^'" configuration due to the cllipticity of 
the image of the lens being very small (and, hence, making it 
more susceptible to having its symmetry axes altered signifi- 
cantly by weak lensing due to foreground galaxies). Naively, 
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Figure 14. Same as Figure 13 except that here the central, ellip- 
tical lens has ehaio = 0.3. Note that the vertical scale differs from 
that of Figure 13. 
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Figure 15. Same as Figure 13 except that here the central, ellip- 
tical lens has ehaio = 0.5. Note that the vertical scale differs from 
that of Figure 13. 



one might hope that the suppression of the observed func- 
tion, {6)/'Y^ (6), could be eliminated simply by choosing 
to compute the mean tangential shear using bright centres 
whose images are highly elliptical. 

In addition, in a search for the signature of anisotropic 
galaxy-galaxy lensing, one might bo tempted to restrict the 
analysis to source galaxies that arc very close to the sym- 
metry axes of the bright centres that are used to calculate 
the mean tangential shear. That is, in all of the analyses 
above, 7'''(^) and 7~(6') were computed using all sources 
whose azimuthal coordinates, ip, placed them within ±45° 
of the symmetry axes of the bright centres. At fixed angular 
separation from an elliptical lens, the maximal difference in 
the shear experienced by two sources will, of course, occur 
when one source is located along the minor axis of the lens 
and the other is located along the major axis of the lens. 
Therefore, one might expect that if one narrowed the anal- 
ysis region from ±45° to, say, ±25° or ±15°, it would be 
easier to detect anisotropic galaxy-galaxy lensing. 

Unfortunately, the situation is not that simple in either 
of these cases. Weak lensing of the bright centres may make 
their resulting images either rounder or more elliptical than 
their intrinsic image shape (i.e., bottom panel of Figure 12). 
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Figure 16. Effect of narrowing the analysis region on the ra- 
tio of mean tangential shears. Here a fiducial elliptical lens with 

(Jv = 150 km scc^^, xt = 100 h^^ kpc, and projected halo ellip- 
ticity ehalo = 0.3 has been adopted. Circles: "Observed" signal. 
Sources have been Icnscd by the central, elliptical lens and all fore- 
ground SIS lenses. The central, elliptical lens has also been lenses 
by all foreground SIS lenses and its observed symmetry axes are 
used for the calculation. Crosses: Sources have been lensed by the 
central, elliptical lens and all foreground SIS lenses. The intrin- 
sic (unlensed) symmetry axes of the central lens are used for the 
calculation. Left: All sources within ±45° of the lens symmetry 
axes arc used in the calculations. Right: Only sources within ±15° 
of the lens symmetry axes arc used in the calculations. Narrow- 
ing the analysis region increases the degree of anisotropy in the 
galaxy-galaxy lensing signal, but it also increases the disparity 
between the observed signal and the one that would be measured 
if the intrinsic symmetry axes of the lenses were known. 



Suppose that one chooses a minimum image cllipticity for 
the bright centres, and that the computation of j^{6) and 
7~(9) is performed using only those bright centres with ob- 
served ellipticity eught > Ccut- Some fraction of the bright 
centres whose intrinsic cllipticity truly exceeds tcut will, by 
weak lensing by foreground galax;ies, have their resulting im- 
ages made rounder than their intrinsic ellipticity. As a result, 
some bright centres with intrinsic ellipticities that axe larger 
than ecut will, in fact, be rejected on because their observed 
(post-lensing) images have ellipticities smaller than Ccut . The 
number of such bright centres that are affected by this will 
vary with the magnitude of the shear that they experience. 

In addition to changing the ellipticity, weak lensing of 
the bright centres may rotate the orientations of their sym- 
metry axes (i.e., top panel of Figure 12). If one simply tries 
to narrow one's analysis region relative to the symmetry itxes 
of the bright centres, a problem will occur if the bright cen- 
tres have been weakly lensed. Any rotation of the syrrmictry 
axes of the bright centres causes the analysis region that one 
truly desires (i.e., the region that brackets the directions of 
the major and minor axes of the projected halo mass) to be 
rotated with respect to the analysis region that one must 
actually use in practice (i.e., the region that brackets the 
directions of the observed major and minor axes of the im- 
age of the bright centre). Therefore, narrowing the analysis 
region may actually increase the discrepancy between the 
observed function, 7^(S)/7^(^), and the function that one 
would measure if the intrinsic symmetry axes of the bright 
centres were known. 

In this section we adopt a fiducial elliptical lens with 
velocity dispersion ay = 150 km sec"'^, truncation radius 



Xt ~ 100 h^^ kpc, and projected halo cllipticity e — 0.3, 
and we construct Monte Carlo simulations that are identi- 
cal to the Monte Carlo simulations in Section 6 (i.e., the 
simulations used to obtain the central panel of Figure 14). 
Shown in Figure 16 is the effect of narrowing the analy- 
sis region when computing the ratio of the mean tangential 
shears. The left panel of Figure 16 shows the observed and 
intrinsic functions, 7 -|-^ 6)/"f~{6), when all sources within 
±45° of the symmetry axes of the central, elliptical lens are 
used for the calculations. The right panel of Figure 16 shows 
the same functions as the left panel, but here only sources 
that are within ±15° of the symmetry axes of the central, 
elliptical lens are used for the calculations. From the right 
panel of Figure 16, it is clear that narrowing the analysis 
region (i.e., using only sources that are very close to the 
symmetry axes) increases the degree of anisotropy in the 
galaxy-galaxy lensing signal. However, narrowing the analy- 
sis region also increases the disparity between the observed 
function, 7^(6)/7~(&), and the function that would be mea- 
sured if the intrinsic symmetry axes of the central, elliptical 
lenses were known. 

Figure 17 shows the effect of rejecting bright centres 
whose images (post-lensing) arc very round. All sources 
within ±45° of the symmetry axes of the central, ellipti- 
cal lens arc used in the calculation. Here circles show the 
observed function, y'^ {0)/'y~ (0), and crosses show the func- 
tion that one would obtain if the intrinsic symmetry axes 
of the central elliptical lens were known. In the case of the 
circles and crosses, no constraint on the ellipticity of the 
lens image is imposed. Triangles in Figure 17 show the ob- 
served function, 7"*" (9) /7~ (9), where the ratio has been com- 
puted using the observed symmetry axes of central, ellip- 
tical galaxies whose images (post-lensing) have ellipticities 
eiight > 0.3. From this figure, then, rejection of lenses with 
image ellipticities eught < 0.3 increases the observed func- 
tion, 7^(0)/7~(6), only slightly. In particular, rejection of 
the lenses with the roundest images docs not allow one to 
recover the function that one would measure if the intrinsic 
symmetry axes of the central, elliptical lenses were known. 



8 SUMMARY AND CONCLUSIONS 

We have investigated the theory of galaxy-galaxy lensing by 
non-spherical dark matter haloes, which should give rise to 
an anisotropy in the tangential shear experienced by dis- 
tant source galaxies. If each distant source is lensed by only 
one foreground elliptical lens, and if the observed symmetry 
axes of the elliptical lens correspond to the intrinsic symme- 
try axes of its projected dark matter halo, one would expect 
the signature of anisotropic galaxy-galaxy lensing to mani- 
fest as 7^ (9) > 7~ {9) over a wide range of angular scales. 
Here 7"*" (9) is the angular dependence of the mean tangential 
shear experienced by sources whose azimuthal coordinates 
place them close to the major axis of the lens, and 7~(S) is 
the angular dependence of the mean tangential shear expe- 
rienced by sources whose azimuthal coordinates place them 
close to the minor axis of the lens. 

Using an observational data set (observed coordinates 
and /-band apparent magnitudes) as a framework for a 
set of Monte Carlo simulations, we have demonstrated that 
the actual signature that one should expect to observe for 
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Figure 17. Effect of rejection of round lenses on the ratio of 
the mean tangential shears. Here all sources within ±45° of the 
symmetry axes of the central, elliptical lens are used. The cen- 
tral, elliptical lens has halo parameters <t„ = 150 km sec~^, 
xt = 100 kpc, and projected ellipticity ehalo = 0-3- Error 
bars are omitted when they are comparable to or smaller than the 
data points. Circles: "Observed" signal. Sources have been lensed 
by the central, elliptical lens and all foreground SIS lenses. The 
central, elliptical lens has also been lenses by all foreground SIS 
lenses and its observed syimnetry axes are used for the calcula- 
tion. No constraint has been placed on the ellipticity of the image 
of the central, elliptical galaxy. Crosses: Sources have been lensed 
by the central, elliptical lens and all foreground SIS lenses. The 
intrinsic (unlensed) symmetry axes of the central lens are used for 
the calculation. No constraint has been placed on the ellipticity of 
the image of the central, elliptical galaxy. Triangles: "Observed" 
signal, computed using only those central, elliptical lenses whose 
observed (post-leiising) image ellipticity is eught > 0.3. Rejection 
of central, elliptical galaxies with very round images docs little 
to affect the discrepancy between the observed function and the 
one that would be obtained if the intrinsic symmetry axes of the 
central, elliptical lenses were known. 



anisotropic galaxy-galaxy lensing is far from the above ide- 
alised case. Because galaxies are broadly distributed in red- 
shift space, it is common for a distant source galaxy located 
at redshift Zs to be Icnsod by another galaxy located at red- 
shift zii < Zs- In turn, this original lens-source pair may 
then be lensed by yet another galaxy (or galaxies) located at 
redshift zi2 < zn. Such instances of "multiple deflections" 
cause the observed signature of anisotropic galaxy-galaxy 
lensing to deviate from the expected signature. The degree 
to which the observed signature of galaxy-galaxy lensing 
deviates from the expected signature is a strong function 
of the characteristic velocity dispersion of the haloes of L* 
galaxies. In the case of low characteristic velocity disper- 
sions, (tJ — 100 km scc^^, the observed ratio of mean tan- 
gential shears, 7^(^)/7~(^), exceeds a value of unity on all 
scales 6 < 100" and is only slightly lower than the function 
one would obtain if the intrinsic symmetry axes of the fore- 
ground galaxies were used to perform the calculation. In the 
case of moderate velocity dispersions, — 150 km sec~^, 
the observed ratio of mean tangential shears shows little to 
no anisotropy on scales > 20". In the case of high ve- 
locity dispersions, a*. — 200 km scc^^, the observed func- 
tion is actually reversed from the expected function (i.e., 
7+(^) < 7~(0)) on scales 20" < 6 < 70", and is consistent 
with no anisotropy on scales 70" < 6 < 120". 



In summary, our simulations show that if one observes 
'y'^{d) ~ 7~(^) in a large galaxy-galaxy lensing data set, the 
observation cannot be simply interpreted as proof that the 
haloes of the lens galaxies are spherically-symmetric. That 
is, although the measured signal appears to be isotropic, 
it is entirely possible that anisotropic galaxy-galaxy lensing 
by non-spherical haloes may have taken place. Further, our 
simulations show that if one observes ^^{0) < 7~(^) in a 
large galaxy-galaxy lensing data set, the observation cannot 
be simply interpreted as proof that mass and light are "anti- 
aligned" in the lens galaxies. That is, although the measured 
signal appears to be reversed from the expected signal, the 
reversal may occur when mass and light are, in fact, perfectly 
aligned within the lens galaxies. 

The primary reason that the observed signature of 
anisotropic galaxy-galaxy lensing differs from the expected 
signature is that the foreground galaxies that are used as 
centres to compute the mean tangential shear have, them- 
selves, been weakly lensed. The expectation that 7^(0) will 
exceed 7~(0) over a wide range of angular scales is based 
upon a picture in which the observed symmetry axes of the 
lenses are identical to the intrinsic symmetry axes of their 
projected dark matter haloes. However, when one computes 
7"'"(&) and 7~ (9) in an observational data set, one cannot di- 
rectly view the intrinsic symmetry axes of the bright, central 
galaxies. Instead, one is forced to use their observed symme- 
try axes and, in general, these will differ from the intrinsic 
symmetry axes. 

Our simulations show that, even in the limit of multiple 
deflections being experienced by the distant source galaxies, 
if one could use the intrinsic symmetry axes of the lenses 
to define the geometry of the problem, one would expect 
to observe 7^(0) > 7~(6). That is, multiple deflections ex- 
perienced by the source galeixies have little effect on the 
intrinsic signature of anisotropic galaxy-galaxy lensing by 
non-spherical haloes. However, weak lensing of the bright, 
central foreground galaxies causes their observed symmetry 
axes (which are used to define the geometry for the calcula- 
tion of 7'''(^) and 7~(^)) to differ from their intrinsic sym- 
metry axes (i.e., the unlensed symmetry axes, which define 
the geometry for the actual lensing of the distant galaxies). 
It is this change in the symmetry axes of the bright, fore- 
ground galaxies that gives rise to the suppression of the ob- 
served function, 7"'"(^)/7~(^), compared to the function that 
would be obtained if the intrinsic symmetry axes were used 
for the calculation. The effects of weak lensing of the bright, 
foreground galaxies on an observation of ^^{0)/^^{9) can- 
not be eliminated simply by rejecting foreground galaxies 
with very small image ellipticities, or by using sources that 
arc particularly close to the observed symmetry axes of the 
foreground galaxies. 

We conclude, therefore, that in order to properly in- 
terpret any observed galaxy-galaxy lensing signal (be it 
isotropic or anisotropic), it is vital that full, multiple- 
defiection Monte Carlo simulations be used. Especially im- 
portant is accounting for the fact that the images of the 
bright, foreground centres are likely to have been weakly 
lensed. If the effects of multiple deflections are not taken 
into account when interpreting an observed galaxy-galaxy 
lensing signal, there is a high probability that incorrect con- 
clusions will be drawn about the nature of the haloes sur- 
rounding the lens galaxies. 
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